Einstein boundary conditions in relation to constraint propagation for the 
initial-boundary value problem of the Einstein equations 
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We show how the use of the normal projection of the Einstein tensor as a set of boundary 
conditions relates to the propagation of the constraints, for two representations of the Einstein 
equations with vanishing shift vector: the ADM formulation, which is ill-posed, and the Einstein- 
Christoffel formulation, which is symmetric hyperbolic. Essentially, the components of the normal 
projection of the Einstein tensor that act as non-trivial boundary conditions are linear combinations 
of the evolution equations with the constraints that are not preserved at the boundary, in both 
cases. In the process, the relationship of the normal projection of the Einstein tensor to the recently 
introduced "constraint-preserving" boundary conditions becomes apparent. 

PACS numbers: 04.25.Dm, 04.20.Ex 



I. INTRODUCTION 



The numerical integration of the initial-boundary value 
problem of the vacuuum Einstein equations is character- 
ized, as is well known [l|, by the presence of four dif- 
ferential equations that must be satisfied by the solution 
of the initial-value problem. These four additional equa- 
tions are normally written in the form of constraints - 
that is: differential equations with no time derivatives of 
the fundamental variables-, and represent the vanishing 
of the projection of the Einstein tensor along the nor- 
mal n a to the time foliation, namely: G a bn b = 0. This 
atypical evolution problem has been a constant source 
of controversy from the onset of numerical implementa- 
tions, because it is commonly found that the numerical 
solution of the initial-boundary value problem with con- 
strained initial data fails to solve the constraints soon 
after the initial time. Earlier numerical simulations pre- 
dominantly opted for a scheme that forced the constraints 
at every time slice (see, for instance, 0). In this scheme, 
the numerical solution of the initial-boundary value prob- 
lem at a given time step is used as a starting guess to 
generate a numerical solution of the constraint equations 
on the same time slice, which is subsequently used as 
data for the initial-boundary value problem for the next 
slice, and so on. In more recent years, this forced con- 
strained scheme has been abandoned in favor of uncon- 
strained simulations that disregard the constraint equa- 
tions, except that the growth of the constraint violations 
is monitored and kept to a minimum by stopping the nu- 
merical simulations before the constraint violations grow 
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too large (see, for instance, 0]). 

The reason for the generically observed drifting of nu- 
merical simulations away from the constrained solution 
has yet to be pinned down and may turn out to be a 
combination of analytic and computational issues. Sev- 
eral analytical factors have been identified that may have 
a role to play in this respect, ranging from ill-posedness 
to asymptotic instability of the propagation of the con- 
straints 0, 0| ■ In the face of the recurrence of un- 
controlled numerical instabilities that prevent numerical 
simulations from maintaining accuracy for a time scale 
commensurate with the emission of strong gravitational 
waves from gravitational collapse, the importance of con- 
straint violations in numerical simulations had widely 
been regarded as secondary until recently. In the 
authors show, with sufficiently accurate numerical exper- 
iments, that the growth of constraint violations plays a 
key role in shortening the running time of numerical sim- 
ulations even when other factors are kept within control. 
Understanding the growth of constraint violations and 
figuring out techniques to keep them in check have be- 
come issues of high priority in the goal to extend the run 
time of numerical simulations to a regime where gravita- 
tional radiation can be extracted. 

What we show in this article is that the projection of 
the Einstein tensor along the normal e a to the bound- 
ary of the simulation, G a f,e b , is equivalent to the con- 
straints that are not identically satisfied at the bound- 
ary by virtue of the initial- value problem alone. This is 
probably a generic fact, irrespective of the details of the 
formulation of the initial-boundary value problem of the 
3+1 Einstein equations. We prove this fact in the case 
of the ADM formulation, which is ill posed but has well- 
posed constraint propagation, and also in the case of the 
Einstein-Christoffcl formulation, which is well posed and 
has well-posed constraint propagation. 
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This fact connects intimately two factors of seem- 
ingly independent effect to numerical simulations; i.e., 
the problem of constraint violations and the problem of 
proper boundary conditions. Intuitively, if boundary val- 
ues travel inwardly at the speed of light (or any finite 
speed), then their effect is confined to the sector of the 
final time slice that lies in the causal future of the bound- 
ary surface, regardless of the extent of the bulk. Thus if 
the bulk is significantly larger than the run time in units 
of the speed of light, the boundary effects will be negli- 
gible in the bulk. As a result, it is commonly assumed 
that boundary effects can normally be kept in check by 
"pushing the boundaries out" and this has largely be- 
come the implicit rule in numerical simulations. (This is 
mathematically true only in the case of strongly hyper- 
bolic formulations, though.) The constraint violations, 
on the other hand, appear to remain in most cases. The 
two problems would thus seem independent of each other 
(in fact, they are itemized separately in |9(). 

On closer inspection, given any formulation of the 
initial-boundary value problem there are constraints that 
are identically satisfied in the simulated region by virtue 
of a proper choice of initial values, and there is a subset of 
constraints that are not satisfied in the entire simulated 
region by virtue of the initial values alone. In a min- 
imallistic sense, it is these "nonpreserved" constraints 
that need to be controlled. Because, as we show here, 
these "nonpreserved" constraints are equivalent (up to 
evolution) to G a i,e b = 0, and because G a be b = 0, in turn, 
are proper and necessary boundary conditions for the 
initial-boundary value problem-as we have shown pre- 
viously 0-, it follows that by treating the boundaries 
properly one may get a handle on keeping the relevant 
constraint violations in check as a byproduct - a reason- 
able deal, by all means. 

The article is organized as follows. The ADM case, 
which is the simplest and most transparent one in spite 
of its modest repertory of useful mathematical attributes, 
is dealt with first in Section^ As an example of strongly 
hyperbolic formulations, the Einstein-Christoffel case is 
developed in Sections IIIII and llVl Both cases are treated 
with the simplifying assumption of vanishing shift vec- 
tor in the terminology of the 3+1 split. The reader may 
assume that the complexity of the argument would in- 
crease sharply in the presence of a non-vanishing shift 
vector, although its general gist should run unchanged. 
This work develops issues that were anticipated in j6|. 



II. ADM FORMULATION WITH VANISHING 
SHIFT 

Throughout the article we assume the following form 
for the metric of spacetime in coordinates x a — (x l , i) in 
terms of the three-metric %j of the slices at fixed value 
oft: 



where a is the lapse function. The Einstein equations 
G a b = for the four-dimensional metric are equivalcntly 
expressed in the ADM form [l|: 



7ij 



-2a K, 



(2a) 



Kij = a (Rij - 2K u K l j + KKij) - D t D,a, (2b) 

with the constraints 
1 



C 



(R - K i:j K ij + K 2 ) = 0, 



(3a) 
(3b) 



Here an overdot denotes a partial derivative with respect 
to the time coordinate (d/dt), indices are raised with the 
inverse metric 7 y , Di is the covariant three-derivative 
consistent with 7^ , is the Ricci curvature tensor of 



R its Ricci scalar, Km is the extrinsic curvature of 



the slice at fixed value of t and K = 7 y . Expressed 
in terms of the Einstein tensor, the constraints © are 
related to specific components in the coordinates (x' l ,t): 



C = -a 2 G u 



C: = - 



OTfij 



G ]t 



(4a) 
(4b) 



where (I4b|l holds only for vanishing shift vector. The 
constraint character (the absence of second derivatives 
with respect to time) is a consequence of the fact that 
the only components of the Einstein tensor that appear 
in |0J have a contravariant index of value t. In geometric 
terms, (@J) are linear combinations of G ao n b = where n b 
is the unit normal to the slices of fixed value of t, and is 
therefore given by n a — g ab n = —ag at — Sf/a. 

Similarly, at any boundary given by a fixed value of 
a spatial coordinate, the normal vector to the boundary 
surface e b can be used to project the Einstein tensor as 
G a be b in order to obtain the components that have no 
second derivatives across the boundary. The vanishing of 
these can then be imposed as conditions on the boundary 
values for the fundamental variables. To fix ideas let's 
choose a boundary surface at a fixed value of x. Thus 
e b = g bx = (0,7") up to scaling. Consequently, G a be b — 
G x a , so any linear combination of the components of the 
Einstein tensor with a contravariant index of value x will 
be suitable. Explicitly we have: 



Gf 

GX 
z 

Gx 



1 



7 " ((ln 7 ), It -7^7^) " KD x a + K x k D k 



K x 1 
-— + KK% + RZ- -D x D v a 
a y v a 



-^l + KK X Z + R x z - —D x D z a 
a a 



ds 2 



-a 2 dt 2 + r )ijdx l dx l 



(1) 



a 

+R1 + - 



1 



{R + K ij Kij +K 2 )+KK X X 



(D j D ja - D x D x a) 



(5a) 
(5b) 
(5c) 

(5d) 
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Here T k i:j = (l/2)7 W (7«j + 7^ - 7^), and the time 
derivative of the components of the extrinsic curvature is 
applied after raising an index, that is: Kj = {^i lk Kkj),t- 
The reader can verify that R x and R x do not involve sec- 
ond derivatives with respect to x of any of the variables 
and that the combination R x — ^R doesn't either. 

At this point the question arises as to whether G a i,e b 
vanish identically on the boundary for any solution of 
the evolution equations © with initial data satisfying 
©• To start with, by inspection one can clearly see that 
Gy and G x are exactly the evolution equations for the 
mixed components K x and K x of the extrinsic curvature 
(namely: Eq. (|2b|l multiplied by 7™ and evaluated at 
j = y and j = z). Thus two of the four boundary equa- 
tions are identically satisfied by the solution of the evo- 
lution equations (irrespective of the initial data). That 
is not the case with G x . If one uses the evolution equa- 
tions for K x and for the trace K that follow from l|2b(l to 
eliminate the time derivatives from G x , one is left with 
the Hamiltonian constraint C. Finally, using the defini- 
tion of the extrinsic curvature (|2a[l into the expression 
for G x it is straightforward to see that Gf = aj xl Ci (this 
is necessarily so in the case of vanishing shift because 
G x = gt a G ax = —a 2 G tx and inverting (j4p|> one has that 
G tx — —a j xt Ci). In summary, if we represent the evo- 
lution equations (|2a|l and (|2b|) in the form £% — ancl 
= respectively by transferring all the terms from 
the right into the left-hand side, we have 

G x t = c£ x + ^^mi - r^diS],) (6a) 

G x = Lfi£* (6b) 

G x = -^Ef z (6c) 

G% = -C + -( 7 fc ^-7^) (6d) 

which can essentially be expressed in the form 

G x ~ aC x (7a) 

G x ~ (7b) 

G x ~ (7c) 

G x ~ -C (7d) 

where the symbol ~ represents equality up to terms pro- 
portional to the evolution equations. Clearly the vanish- 
ing of G x and G x at the boundary is intimately connected 
with the propagation of C and C x towards the boundary. 
On the other hand, C y and C z are not related to the 
boundary equations. 

To have a sense for why this should be so, one may 
now look at the auxiliary system of propagation equa- 
tions for the constraints as implied from J5J. Taking a 
time derivative of the constraints and using the evolution 
equations we have 

C = ad l C l + ... (8a) 
Ci = ad t C + ... (8b) 



where . . . denote undifferentiated terms. This is a well- 
posed system of equations in the sense that it is strongly 
hyperbolic 0. The characteristic speeds are and ±<x 
With respect to the unit vector = j ix j \J^ XX , normal 
to the boundary, the characteristic fields that travel with 
zero speed are C v and C z . This is an indication that 
these constraints could be interpreted as "static" , and 
that their vanishing is preserved by the evolution at all 
times after the initial slice, that is: they will vanish at 
the boundary by virtue of the evolution equations and 
the initial values (of course, only so long as the other 
constraints also vanish, due to the coupling in the un- 
differentiated terms in the constraint propagation equa- 
tions). It makes sense, thus, that they do not relate at 
all with the boundary equations G x = and that two of 
the boundary equations are redundant (i.e., trivial). 

The characteristic fields that travel with non-zero 
speeds ±a are, respectively, C ± = C ± C x /^/j xx . The 
characteristic field C+C x / y/j xx is outgoing at the bound- 
ary, which means that its value is propagated from the 
initial slice to the boundary. If this constraint combina- 
tion vanishes initially, then it should also vanish at the 
boundary. On the other hand, the characteristic field 
C — C x I \J~j xx is incoming at the boundary, which means 
that its value at the boundary is unrelated to the initial 
values prescribed in the interior. Thus, this combination 
of constraints is not vanishing at the boundary by virtue 
of the initial data and the evolution. This is a "non- 
preserved" constraint. This constraint must be enforced 
at the boundary by means of appropriate boundary con- 
ditions if the solution of the ADM equations @ at the 
final time slice is to satisfy all four of the constraints. 
Since the two boundary equations G x = and G x — 
are directly related to the two characteristic constraints 
with non-vanishing speeds, imposing either one or any 
linear combination of them except G x + (a/ \/j xx )G x = 
on the boundary is equivalent to imposing a linear rela- 
tionship C~ = BC + (whith a constant B) as a boundary 
condition on the system of equations for the constraints. 
These are valid boundary conditions that preserve the 
well-posedness of the evolution of the constraints Q- 

Thus, of the two equations G x — and G x = 0, one 
of them is necessary to avoid constraint violations. The 
remaining one, which we can represent by the linear com- 
bination G x + (a/ \/j xx )G x — 0, is equivalent to a con- 
dition on the outgoing constraint C + . Therefore, it is re- 
dundant to the system of evolution of the constraints. At 
this point, however, it is not clear whether this remain- 
ing equation would be redundant as a boundary condi- 
tion for the ADM evolution. The reason is that since the 
ADM equations are not strongly hyperbolic, the number 
of boundary conditions it requires is not known. 

We are led to infer that some of the components of the 
projection of the Einstein tensor normally to the bound- 
ary are in direct correspondence with those constraints 
that are not automatically preserved by the evolution 
equations and the initial values. In this particular case 
of the ADM formulation, because the evolution equations 
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thenselves are not strongly hyperbolic, perhaps nothing 
else can be said about the role of the equations G a be b = 
as boundary conditions, except that failure to impose the 
nontrivial one leads to constraint violations with guaran- 
teed certainty. 

Yet the case of the ADM equations has been sufficiently 
direct and transparent to provide us with a guide for the 
much more complicated (but also much stronger) case of 
formulations of the Einstein equations that are strongly 
hyperbolic. In the next Section, a particular formulation 
is chosen for its simplicity and its close relationship to 
the ADM case. 



III. BOUNDARY CONDITIONS FOR THE EC 
FORMULATION WITH VANISHING SHIFT 

Most hyperbolic formulations of the Einstein equations 
require the introduction of additional variables in order 
to reduce the differential order in the spacelike coordi- 
nates from second to first. The additional variables are 
always a complete set of linearly independent combina- 
tions of the space-derivatives of 7y. This inevitably in- 
troduces an ambiguity: whereas indices labeling second- 
order derivatives are automatically symmetric (which af- 



7y 

Kij + aj kl dif ki] 



f kij + ad k Kij 



with the constraints 

C = -\l ij l kl {2{d k f i j l -d i f jkl )+K ik K jl -K ij K kl 

+Y nn [fikm(5fjln — Gfljn) + ^fiklfjmn 
+fijk(Sfmln ~ 20f lmn }} = (11a) 
Ci = —7 kl {Y nn [Kik(3flmn — 2fmnl)—Kkmfiln] 

+diK k i-d k K u }=0 (lib) 

Ckij = 2fkij - ^ lm {flm{il])k - lk{ifj)lm) - dk^ij = 

(lie) 

to be imposed only on the initial data. Here C^ k rep- 
resent the definition of the additional 18 first-order vari- 
ables fkij needed to reduce the ADM equations to full 



fords plenty of convenient cancellations) , the same terms 
written in terms of the new variables are not manifestly 
symmetric, the symmetry becoming a commodity that 
may or not be "turned on" by imposing yet additional 
constraints. For example, in second order one has, un- 
ambiguously, "fij^ki — lij.ik = 0. However, if one defines 
dijk = lij.k, one is free to write the same expression ei- 
ther as or as dij k j — diji t k- Undoubtedly, it makes a 
big difference to a numerical code to write it one way or 
the other, since they are not the same when it comes to 
the solution-generating process. This ambiguity is usu- 
ally taken advantage of in order to write the Einstein 
equations in manifestly well-posed form. The cost is the 
loss of some transparency and the addition of constraints 
and, as we will see shortly, boundary conditions. 

Consider the Einstein-Christoffel [EC] formulation of 
the 3+1 equations §. By defining first-order variables 
as the following 18 linearly independent combinations of 
the space-derivatives of the metric: 

fkij = ^(ij)k + lkil lm ^{lj\m + lkjl lm ^[li\m, (9) 

and choosing the lapse function as a = Q^/7 with Q as- 
sumed arbitrarily prescribed a priori, the 3+1 equations 
can be put in the following equivalent form 0: 



(10a) 



(10b) 



(10c) 

I 

first-order form (they represent © inverted for 7y,fc in 
terms of fkij)- No other constraints are needed to pick 
a solution of the Einstein equations out of the larger set 
of solutions of the EC equations. A point that deserves 
to be made here is that l|10b|l is an exact transcription 
of the ADM evolution equation J2EJ, with no mixing of 
constraints into it. In other words: Eq. llObfl reduces 
exactly to Eq. (|2bll if fkij is substituted back in terms of 
jij (assuming the same lapse condition). 

We need to translate the boundary equations JSJ in 
terms of the EC variables. Because of the ambiguity 
in writing second derivatives 7^^ = 7^ = (jij.ki + 
7y./fc)/2 in terms of first derivatives of fkij there are, 
in principle, many different ways to write — in 



= -2aKij 

= a{j kl (K k lKij - 2K kt Kij) + 7 ^ mn {^fkmif[lri[j + ^fkm[nft]ij ~ fikrnfjln 
+&f(ij)kf[ln]m + 

m{ifj)ln &fklifmnj ~\~ 20 fkl(ifj)mn ^^fiklfjmn) 

-d^ InQ-di In Qdj In Q + 2 7ijl kl 1 mn (f kmn d l In Q - f kml d n In Q) 
+l kl [{2f(i j)k ~ fkn)dl In Q + ifkliidj) In Q - 3(f m d 3 In Q + f jkl d t In Q)}} 

= a{j mn [4K k ^fj) mn - ^fmn{i K i)k + Kij(2fmnk ~ 3fkmn)] 

+ 2y nn ~f Pq [K mp ("fk(ifj)qn - 2f qn ( l Jj)k) + J k (iKj)m(8fnpq ~ Gfpqn) 

+K mn (4f pq( aj) k - 5"f k (ifj) pq )} - Kijd k \nQ 
+2j mn (K m(l j j)k d n \nQ - K mn -i k{i d 3) InQ)} 
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terms of fkij- Notwithstanding, among all these differ- 
ent possibilities, the admissible boundary equations are 
determined by the requirement that no x— derivatives of 
any variables may occur. In the following, we write the 
simplest expression of the boundary equations JSJ up to 
undifferentiated terms. 

We start with Eq. IJ5EJ. The terms with second deriva- 
tives arise from the combination R x — j xk a, y k /a where 
a = Qy/j with arbitrary Q. Up to undifferentiated terms 
in the fundamental variables this combination is 



R. 



x 7 a iyk _ 1 x 



= —-^7 Xm l (lym,M — Jyl,km ~ lkm,yl) 



n ,xm n ,kl, i 
"7 7 lkl,ym + 



(12) 



which can equivalently be represented as 



x "j x a, y k 1 xm kl 

Ry ~ 



n = 7 Xfym,kl — lyl,km) 



+ o7 7 Uylkm,l 



1 

2 

-7 7 OyJkl,m + 



(13) 



where it is clear that in the last two terms jkmj and ^ki,m 
can be substituted now directly in terms of fkij via 

lijM = 2fkij - 4:l lm (flm{il])k - 7k(ifj)lm) (14) 

yielding 



Using all these in (|13fl we finally have 
R x v - = d z (/V - /V) 

+dy (f X y V ~ f V y X + 3/ V - 2/ V) 

+ ... (21) 

The principal terms in R x ~j xk a, z k /a can, evidently, be 
handled in the same manner, so we have the following ex- 
pressions for (|5b|l and l(5cjl explicit up to undifferentiated 
terms in the fundamental variables: 

Gy = -^+M/V-/V) 

+9y (f X y V - f V y X ~ 3/ X k" + 2/ V) +■■■ (22) 

G x = -^+d y (f x z y -f y z x ) 
a 

+d z (/V - /*/ - 3/V + 2/V) + ... (23) 

The calculation of the representation of G x in terms 
of fkij is much longer, but equally straightforward. In 
the following only the main steps are indicated. We start 
with the explicit expression of the Ricci components in 
terms of the metric: 

ryx 1 r> ^ {„,xm„kl „,xl„,krn\„, 

1 



77 (7 7 -7 7 llymM 



^l xm l kl d ylkm ,i = fy(-3/V + 4/V) + ... (15) 



l xm l kl d ylkl , m = d y { l/V • 6/V) • ... (16) 

Thus the last two terms clearly have no x— derivatives. 
By a relabeling of the dummy indices the remaining term 
is equivalent to 



7 7 \lym,kl~lyl,km) = 17 7 _ 7 7 )lym,kl (17) 

One can see by inspection that the contribution of k = x 
is identically vanishing, so we have, equivalently: 

(„,xm^kl „xl „,kvn\„ / (~,xm^ s yl „,xl „,ym\ , 

(7 7 -7 7 hymM = (7 7 -7 7 )7am,j,i 

+ (7 7 -77 nym,zl 

= Oy [(7 7" - 7 7" )7ym,iJ 

+9, [( 7 Km 7 z ' - 7^7^)7^,;] 
+ ... (18) 

Substituting 7 ym ,i in terms of fkij via (|14|) we have 



-j(7 7 -7 7 )7zm,fei + ■ ■ • ( 24 ) 

Expanding appropriately some of the contractions indi- 
cated (with the guidance that no second a;— derivatives 
may remain) leads to a large number of cancellations, 
yielding the equivalent form: 

K-\R = -\{l Zm l yl -l zl l ym )ly m ,zl 

-\(l ym l zl -l yL l zm h zm ., y i + ...(25) 

where it is obvious that d z and d y can be pulled out of 
terms that are expressible in terms of fkij ■ We have 



ii zm i yl - 7 V m hw = /V-/V 



{i ym i zl -i yl i zm )iz m ,i = /V-/V 



+f z k k ~ f k k z + ■ • ■ (26) 
+/V - f k k y + ■■■ (27) 



{l xm l yl -l xl l ym )ly m ,i = /V /V 

-/V + r fe" + ■ • ■ (19) 

(fY-fVlv,! = /V-/V + -- (2°) 



Consequently we have 



1 



^ - 2^ = 0»(/V - /v + /v - /V) 

+d y (f y z z - f zy + f k k v - Fk k ) + ■ ■ -(28) 
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The terms on second derivatives of the lapse in (|5dl 
must be represented in terms of fkij as well. They are: 

-(D k D k a-D x D x a) = jV m ^y + + 
a a a 



X7 (7 lkl,my+l lkl,mz) 



2/ R fe 



(3/V - 2/V) 



(29) 



With and the boundary equation J5dj reads 



Gx 



K - K% 
+d y (f V y V 



2/V - 
■2/V 



J £C ' 

fxy 
J x 



Zf Z y V 
- 3/V 



2/V) 



2/V) 



(30) 



Finally, it is quite straightforward to translate the 
boundary equation Ij5a(l in terms of fkij '• 



ih 3 

J k 



(31) 



where /V = d t (j km f„ 

We thus have the four equations (|3*T|) . (j3UJ|, an d 

1|23|1 to be considered as potential boundary conditions 
for the EC equations. Which ones are identically satisfied 
by virtue of the EC equations and the initial values, and 
which ones remain as boundary conditions is determined 
by the characteristic fields of the EC equations. 

The EC equations are symmetric hyperbolic (and thus 
strongly hyperbolic as well [?]])■ With respect to the unit 
vector £ l = j xl / yj^ xx which is normal to the boundary 
x ~ xo for the region x < xo there are 18 "static" char- 
acteristic fields (the six 7^ , the six f y ' j and the six /'V) 
and 12 characteristic fields traveling at the speed of light, 
of which six are incoming: 



-U? = K 



fx 
J i 



and six are outgoing: 



■w 



(32) 



(33) 



With regards to boundary values, those of the static and 
outgoing fields are completely determined by the initial 
values. Because the time derivatives of the incoming 
fields ~Uy 7 ~U* and the combination ~Ujj — ~U Z do 
not occur in any of the four boundary equations <|3"T|) . 
()30JI . (|22|l or (|23(l . the values of these three incoming 
fields on the boundary are completely arbitrary. How- 
ever, the time derivatives of ~Uy and occur in i|22|) 
and respectively, so the boundary values of ~Uy 

and ~U~z are determined by (1221 and (|23|l in terms of 



outgoing and static characteristic fields, as well as initial 
values. Finally, the time derivatives of the combination 
~U-j + ~U Z appear in (|3"T|l and in (j20J), as well as those 
of the corresponding ougoing counterpart + U-j + + U Z . 
This means that the two equations (|31|l and (|30|l are 
equivalent to one boundary prescription for an incoming 
field {~U-j + ~U Z ) and one condition on an outgoing field 
( + U y + + U Z ) which must be identically satisfied. 

Thus, in the EC formulation of the Einstein equations, 
three of the four equations G a be b = are non trivial, as 
opposed to one in the ADM case. 



IV. RELATION TO CONSTRAINT 
PROPAGATION IN THE EC FORMULATION 
WITH VANISHING SHIFT 

The question that arises in the EC case is how the three 
non-trivial boundary conditions in the set G a b& h = re- 
late to the constraints. In analogy with the ADM case of 
Section[n] we anticipate that in the EC case there will be 
(at most) three constraints that will be "non-preserved" 
at the boundary, and that they will be related to the three 
non-trivial components of G a be b = by linear combina- 
tions with the evolution equations. To prove this would 
be far from trivial, in principle, but can indeed be done 
explicitly by using the ADM case as a guide, as we show 
in the following. 

We start by looking at the auxiliary system of propa- 
gation of the constraints. The reader should note, in the 
first place, that the addition of the 18 first-order con- 
straints Ckij automatically raises the differential order of 
the system of propagation equations for the whole set of 
22 constraints to second order: the evolution equations 
for the set (C,Ci,Ckij) implied by the EC equations in- 
volve second space-derivatives of Ckij ■ We have explicitly: 



C 

c\ = 

Ckij — 



ad l C 
1 



d k {d i C kl l -d l C ki l + d k C li l -d k C il l ) 
-ad,C + ... 



(34a) 

) 

(34b) 
(34c) 



where . . . denote undifferentiated terms. The system of 
auxiliary evolution equations l|34|l for the 22 constraints 
is not manifestly well posed because of the presence of 
second-derivatives in the right-hand side (in fact, some 
times this indicates ill-posedness 0])- However, it can 
be reduced to first differential order in the usual man- 
ner, by adding new "constraint" variables that are space- 
derivatives of the constraints. In this case we can do with 



Clkij = 7J i^lCkij — dkClij) 



(35) 



which enlarges the system l|34|) to 40 variables in all, and 
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casts it in the following form: 

C = ad l C l + ... (36a) 
Ci = -a (diC + d k C lkl l + d k C ka l ) + . . . (36b) 
Cki 3 = ... (36c) 
Cikij = ailkidiCj +j k] diC i - -fudkCj ~ jijd k Ci) + . . . 

(36d) 

We digress slightly now from the main point to point out 
that, if one chooses to do so, the new "constraint" vari- 
ables Cikij can be expressed in terms of the fundamental 
variables of the EC system, in which case they read: 

Cikij = dl(fkij — 2i nm (fnm(ilj)k — Ik(ifj)nm)) 

-dk(flij - ^I nm {fnm(ilj)l - ll(ifj)nm)), (37) 

and turn out to be identical to the "integrability condi- 
tions" = l/2(di n /ij i k — dk7ij,i). This is usually the way 
that such "constraint quantities" necessary to reduce the 
propagation of the constraints to first differential order 
are introduced by most authors, starting with Stewart in 
[Tl|. They are not to be interpreted as additional con- 
straints to impose on the initial data, though. In fact, 
Cikij = holds identically for any initial data satisfying 
Ckij = 0. They are not to be used as substitutes for 
Ckij = either, because the vanishing of Ckij does not 
follow from the vanishing of Cikij ■ Indeed, if Cikij = 
is satisfied by 7^ and fkij , then there exists a field hij 
such that 2 - 4rf lm {fi m {i1j)k ~ lk{ifj)lm) = d k hij, 
but it does not follow that hij — 7^ . The introduction 
of the additional "constraints" Cikij — is mostly a for- 
mal procedure to study the evolution of the constraints 
as functions of the point -not through the fundamental 
variables. In fact, one does not need to include all of Cikij 
as additional variables in order to reduce the propagation 
of the constraints to first differential order, since one can 
see that only the contractions Ciki 1 and Cku 1 appear in 
the equations, and these are 12 variables, so at least six 
of the 18 new variables Cikij are entirely redundant. 

The immediate benefit of introducing the first-order 
constraint variables is that the system (|36J) is well posed 
in the sense that it is strongly hyperbolic with charac- 
teristic speeds of (multiplicity 34), +a (multiplicity 3) 
and — ol (multiplicity 3). This means that, with respect 
to the (outer) boundary at x = Xq, three characteristic 
constraint variables are incoming, three are outgoing and 
all the others are "static" . The six non-static character- 
istic constraints, which we denote by ± 2i (with + for 
outgoing and — for incoming), are explicitly: 

±Z X = C X ±^L=(C + C kx xk +C x xk k ) (38a) 

± Zy = Cy±-j== (C kX y k Hh C* yf) (38b) 

± Z Z = C z ± -±= (C kx zk + C\ k k ) (38c) 

Because three of the characteristic constraints are incom- 
ing at the boundary, even if they are set to zero initially 



they will not be vanishing at the boundary by virtue of 
the evolution equations. In fact, the three incoming con- 
straints must be prescribed at the boundary, either arbi- 
trarily, or as functions of the outgoing constraints. That 
is: the problem of the propagation of the constraints, 
Eq. H36fl . requires three boundary conditions. Thus we 
have proven the first part of our claim: there are three 
constraints that are "non-preserved" at the boundary, in 
the same number as nontrivial boundary conditions for 
the EC equations, as anticipated. 

The second part of the claim is to prove that the incom- 
ing constraints are related to the equations G x = by 
terms that are proportional to the evolution equations. 
The size of the problem as regards the number of vari- 
ables makes this practically impossible unless one had 
a good guess as to what the linear combinations ought 
to be. Guided by the ADM case, we may expect that, 
through the evolution equations, G x will be related to C 
and G x will be related to C x , except for terms involving 
the new constraints Ckij or Cikij ■ We also expect that G x 
and G x will only be related to the new constraints, but 
not the hamiltonian nor vector constraints. Our ansatz 



is explicitly as follows: 

G x ~ aC x (39a) 

G x ~ C kx yk +C x yk k (39b) 

G x ~ C kx zk +C x zk k (39c) 

G x ~ C + C kx xk +C x xk k (39d) 



which is equivalent to © except for terms whose occur- 
rence is suggested by the combinations of constraints that 
can be expressed directly in terms of linear combinations 
of the characteristic constraints ± Z;. 

The ansatz may be verified (or disproven) simply by 
comparing the right hand sides of equations (|31|l . J22J, 
(|23|l and (|3U|) with (|39l) term by term under the assump- 
tion that every occurrence of a time derivative in 13 1L 
(|22p. (|23|l or l|3U|l must be substituted in terms of space- 
derivatives by means of the evolution equations, that is: 

Kf = <u) l ,f k / • ... (40a) 
/V = -ad k Ki + ... (40b) 

In particular, using l|40a(l to eliminate K x , for 1221) we 

have 

G x ~ dify 1 + d y (2f k k x - 3/V'i + • ■ • (41) 
On the other hand, by definition from l|37(l we have 

C^yk +C X y k k = dlf X y l +dy(2f k k X ~ 3/^) + •-. (42) 

Thus 1)39 b|l is verified. By the same argument substitut- 
ing y with z, (|39c|) is similarly verified. 

Next, using (jllajl to eliminate K — K%{= K^j + K z ), 
for l|30|) we obtain: 

(■:;. - -d*(/y + /v) 

+d y (2/ V - /V + 2/V - 2/V) 

+d z (2/V - /V + 2/V - 2/V) + • • • (43) 
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On the other hand, directly by the definition 137|) and 
the expression (|llaf) for the hamiltonian constraint we 
have 

r> , nkx , nx k 
U T U x k T" a;fc 

= -^(/V + /V) 
+d y (2/v - /v + 2/v - 2/V) 

+9 Z (2/%^ - /V + 2/V - 2/»/) + . . . (44) 

Thus also l|39d(l is verified. 

Finally, using (|40b|) to eliminate the time derivatives 
of /"V m terms of space derivatives of if?, Eq. reads 

G x t ~ - a (d x K -d k K£) + ... (45) 

which is manifestly equal to aC x , thus verifying (|39a|) . 

If we represent the evolution equations IjlObjl and (|10c|) 
in the form £y = and ffc^ = 0, respectively, by simply 
transferring all the terms in the right-hand side to the 
left, what we have shown is that the following relation- 
ships between the projection of the Einstein tensor, the 
constraints and the evolution equations hold in the EC 



formulation: 

G x t = aC x +£ x k k -£ k k x (46a) 

G x = C kx yk + C V - (46b) 

G x = C kx zk + C x zk k - a~ x ei (46c) 



G x = C+C kx xk +C x xk k + a- 1 £^~a- 1 £ x (46d) 

It may be objected that the proof is not complete be- 
cause the undifferentiated terms of the equations have 
not been shown explicitly to be the same on both sides. 
We think that it is quite clear that a complete proof in 
that sense may not be feasible, but yet the consistency of 
the principal terms and the geometrical foundation of all 
terms on both sides, taken together, give a very strong 
indication that the equality will hold term by term. 

The relationship to the characteristic constraints is 
found by "inverting" (|38[) in order to have the fundamen- 
tal constraints in terms of the characteristic constraints: 

Ci = ^(+Z i + -Z i ) (47a) 
C kx yk +C x yk k = ^L(+ Zy --Z y ) (47b) 
C kx zk +C x zk k = ^1(+Z Z --Z Z ) (47c) 

C + C kx xk +C x xk k = ^1(+Z X --Z X ) (47d) 

Thus, up to terms proportional to the evolution equa- 
tions, imposing G x — is equivalent to imposing 



1 x \+Z l + -Z. i ) = (48a) 

+ Z y --Z y = (48b) 

+ Z Z --Z Z = (48c) 

+Z X --Z X = (48d) 



From our discussion in Section II I II it follows that (|48b[) 
and H48cf> need to be imposed, and then either (|48all or 
(I48d|l or, in fact, any linear combination of them, as the 
third boundary condition for the EC equations. Remark- 
ably, this is entirely consistent with the constraint prop- 
agation problem l|H6jl. since l|48b|l . I|48c|) plus one linear 
combination of (|48a|l and i|48d|) turn out to be an admis- 
sible complete set of boundary conditions for the three 
incoming constraints in terms of the three outgoing con- 
straints. By admissible boundary conditions we mean, 
for instance, maximally dissipative boundary conditions, 
which, as is known 7], for homogeneous problems like 
this one, essentially take the form of ~U{ = L^Uj with 
a rather general (bounded) matrix L\ . 

In fact, the reader will have no difficulty to recog- 
nize in l|48a|) , l|48b|l and (|48cl) the "constraint-preserving" 
boundary conditions of the Neumann type of |12| (with 
rj = 4), where boundary conditions arising from con- 
straint propagation are discussed in a linearized setting 
(and furthermore restricting r\ to the interval < rj < 2.) 
This means that the set of three boundary equations 

G x = 0, 
G x y - 0, 
G x = 0, 

for the EC formulation -with G X ,G X ,G Z given by J5a}, 
(|5bp and (|5c|l written in terms of the variables fkij in 
such a way that the principal terms occur as in 13 H . 
(|22|l and 123|) respectively- constitute the exact (non- 
linear) form of the so called constraint-preserving bound- 
ary conditions of the Neumann type that one would write 
for the case of the Einstein-Christoffel formulation by 
following the prescription in This result general- 

izes to three dimensions the prediction that we advanced 
in |13| . where the constraint preserving boundary con- 
ditions were found to be identical to the projections of 
the Einstein tensor normal to the boundary for the EC 
equations with the restriction of spherical symmetry. 

What needs to be made clear, however, is the ques- 
tion of why should the equations G x — imposed on the 
boundary enforce the constraints in the interior, given 
that they are only equivalent to the constraints at places 
where the evolution equations are satisfied. As a mat- 
ter of fact, the evolution equations are not imposed on 
the boundary (that's why one needs boundary prescrip- 
tions). Therefore it is not true that the constraints are 
being enforced on the boundary. Still, they are enforced 
in the interior, which is the goal. The argument is the 
following. The vanishing of G x on the boundary is, 
in a sense, "carried" by the incoming fields to the inte- 
rior, where the evolution equations are actually satis- 
fied. In the interior, then, one has both G x = and 
£ij — £ij k = 0. Therefore, by i|4l)|) . the constraints re- 
lated to G x are enforced in the interior. Thus, by impos- 
ing G x = on the boundary, one enforces the constraints 
that would otherwise be violated outside of the domain 
of dependence of the initial surface. 
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V. STABILITY CONSIDERATIONS 

The fact that three linearly independent combinations 
of G x — except Gf + -f xl Gf — yield necessary and 
sufficient boundary conditions on the incoming funda- 
mental fields of the EC equations has importance in its 
own right: they ensure the uniqueness of solutions to the 
initial-boundary-value problem. Additionally, the fact 
that such boundary conditions are related to the non- 
static characteristic constraints by linear combinations 
with the evolution equations is entirely equivalent to a 
procedure of "trading" transverse space derivatives by 
time derivatives using the evolution equations, which is 
the procedure that has been used in |12] to obtain what 
the authors refer to as constraint-preserving boundary 
conditions. Which of the two meanings (linear combina- 
tions with evolution versus trading derivatives) to use is 
only a matter of taste. We (as the authors of 0] a P~ 
parently do) interpret this fact as an indication that the 
boundary conditions of the evolution problem are equiva- 
lent to imposing vanishing values of the constraints on the 
boundary (up to linear combinations of the constraints 
among themselves). 

The reader with an interest in numerical applications 
may be concerned with the question of whether any or 
all sets of Einstein boundary conditions (by which we 
mean three linearly independent combinations of G x a = 
except Gf + "f xl G x = 0) preserve the well-posedness of 
the initial value problem represented by the EC equa- 
tions. More precisely, with the question of which sets of 
Einstein boundary conditions ensure that the solution at 
the final time depends continuously on the initial data. 
This question has not been answered in the previous sec- 
tions, as it lies beyond the scope of the present work, 
requiring a calculation by means of either energy meth- 
ods (in the best case) or Laplace transform methods (in 
the most likely case), as described in 0. 

However, preliminary results of relevance to this ques- 
tion already exist in the recent literature, and it is ap- 
propriate to mention them and their immediate con- 
sequences. In jf4| the authors undertake a Laplace- 
transform type of study of the well-posedness of sev- 
eral boundary prescriptions including some choices of the 
Einstein boundary conditions for the standard EC sys- 
tem, as in the pres ent article (as opposed to "general- 
ized" EC as in 12]). In Section IV of Q the authors 
conclude that, in the linearization around Minkowski 
spacetime, the choice of three boundary conditions as 
G x t = G xy = G xz = is in fact well posed for the case 
of standard EC as in the present article. Since the lin- 
earization of Gf = Gy = G x = around Minkowski 
space coincides with G x t — G xy = G xz = 0, this in fact 
proves that, in the linearization around flat space, the 
three Einstein boundary conditions Gf — G x = G x — 
constitute a well-posed set of boundary conditions for the 
(standard) EC formulation. 

Moreover, by the present paper, these three Ein- 
stein boundary conditions represent the well-posed Neu- 



mann boundary conditions of the "constraint-preserving" 
scheme of |l2| if one were to extend those authors' argu- 
ment and lexicon to the case rj = 4. In the constraint- 
preserving scheme of 0] the authors did not allow the 
parameter of the generalized EC system to take the value 
rj = 4 required for the standard EC system, for the rea- 
son that the energy method of works well only for 
symmetric hyperbolic formulations with symmetric hy- 
perbolic constraint pro pag ation, which is not the case for 
-q > 2. As a result, in [13] the authors prove that rj < 2 
(necessary for symmetric hyperbolic constraint propa- 
gation) is sufficient for well-posed constraint-preserving 
boundary conditions, but they do not prove that rj < 2 
is necessary. On our part, by combining the results 
of ^3] (in which the Laplace-transform method is used 
rather than energy methods) with our current article 
as indicated, we are indeed demonstrating that r\ < 2 
is, in fact, not necessary in order to have well-posed 
constraint-preserving boundary conditions of the Neu- 
mann type (this fact is actually indicated implicitly - 
but unambiguously- in ^3])- This also disproves, by 
counterexample, the erroneous belief that the propaga- 
tion of the constraints needs to be symmetric hyper- 
bolic (not just strongly hyperbolic) in order for such well- 
posed constraint preserving boundary conditions to exist. 
But more interestingly, we are providing the principal 
terms of the exact nonlinear boundary conditions, the 
non-principal ones being precisely indicated by the ex- 
act expression of the components of the Einstein tensor, 
namely by full substitution of the fundamental variables 
in Eqs. ©. 

It also needs to be p ointed out to the readers that 
in the same paper 1 1 -11 ] the authors show that impos- 
ing the linear combination G xx — Gtt = in addition 
to G xy — G xz — leads to an ill-posed problem in the 
linearization around Minkowski space. This is an indi- 
cation that Gf - j xl Gf = used with G x = G x = 
may lead to an ill-posed initial-boundary-value problem 
in the nonlinear case as well, though an actual proof for 
the nonlinear case is still lacking. As terminolog y g oes, 
readers should be made aware that the authors of [Tj] in- 
appropriately associate the concept of Einstein boundary 
conditions exclusively with this set, without acknowledg- 
ing that, in the way we have used it, the term refers to all 
linear combinations of G x = 0, including the well-posed 
ones, and including those that, as we show, are equiva- 
lent to the "constraint-preserving" boundary conditions 
of [2 in the linearization. 



VI. CONCLUDING REMARKS 

What has been shown is that the vanishing of the pro- 
jection of the Einstein tensor normal to the boundary 
acts as boundary conditions for the EC equations that 
automatically enforce the constraints that are otherwise 
not propagated, and that the set Gf — G x — G x z = 0, in 
particular, leads to a well-posed initial-boundary-value 
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problem ideally suited for numerical evolution. 

The reader can see that the concept of Einstein bound- 
ary conditions can be applied to any formulation of the 
initial-boundary value problem of the Einstein equations. 
The choice of the EC formulation in this work is merely 
done for convenience as the clearest illustration. We 
think it is reasonable to postulate that in the case of any 
strongly hyperbolic formulation of the Einstein equations 
the Einstein boundary conditions will include a set that 
will be equivalent to constraint-preserving boundary con- 
ditions of the Neumann type. But most numerical simu- 
lations today are currently done using formulations that 
are not strongly hyperbolic, in which case the question 
of well-posedness becomes irrelevant, yet the problem of 
identifying useful and consistent boundary conditions re- 
mains. For such formulations, the Einstein boundary 
conditions eliminate some guesswork on boundary values 
and may help control the constraint violations. 

Given the explicit expressions (JSJ , the Einstein bound- 
ary conditions G a be b — (up to linear combinations as 
indicated throughout) for any 3+1 formulation of the 
Einstein equations are found by simply expressing JSJ in 
the chosen representation of fundamental variables. The 



only subtlety that must be taken care of is that the final 
form of G a i,e b — in the chosen variables must contain 
no derivatives across the boundary of the same order as 
the evolution equations. That such a form exists is guar- 
anteed by the Bianchi identities. In all cases, the proce- 
dure to obtain the correct form of the Einstein boundary 
conditions in the chosen variables parallels what is done 
in the current work for the case of the EC formulation. 
In particular, with perhaps very little effort our method 
can almost certainly be used to show (or disprove) our 
claim that the constraint-preserving boundary conditions 
of the Neumann type for the generalized EC systems as 
in |l2j| are, indeed, Einstein boundary conditions as well. 
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